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conduction  problem  is  presented.  A  greater  benefit  is  obtained  by  increasing  the 
number  of  nodes  rather  than  decreasing  the  integration  time  interval.  A  procedure 
is  given  for  matching  the  predicted  and  measured  initial  values  of  the  deep  body 
temperature  and  the  metabolic  rate  which  allows  a  more  accurate  determination  < 
set-point  values  for  thermoregulation.  To  circumvent  the  acute  sensitivity  to  th 
6kin-water  temperature  difference  using  the  conventional  prediction  of  convective  hea 
loss,  use  of  a  heat  balance  during  steady-state  of  the  skin  temperature  is  outlined. 
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Nude  immersion  in  water  colder  than  the  body  temperature  represents  an 
acute  exposure  to  cold  since  convective  heat  loss  is  many  times  greater  than  in  air. 
Temperature  gradients  become  large  and  physiological  responses  are  dramatic. 
These  factors  contribute  to  the  complexity  of  mathematically  modelling  human 
thermoregulation  for  cold  water  immersion.  Some  recent  attempts,  though,  have 
been  made  with  reasonable  success.  These  include  multi-segment  models  that  use 
physical  representations  of  the  human  body,  principles  of  heat  conduction  and 
feedback  control  mechanisms  (Miller  and  Seagrave  1974;  Montgomery  1974;  Wissler 
1985).  In  an  effort  to  improve  the  predictive  capability  of  such  models,  three 
specific  problems  will  be  addressed:  1)  finite-difference  solution  of  the  bioheat 
equation;  2)  differences  between  predicted  and  measured  initial  conditions,  and  3) 
prediction  of  the  convective  heat  loss. 

The  basic  multi-segment  thermoregulatory  model  can  be  separated  into  two 
distinct  systems,  the  controlled  system  and  the  controlling  system.  The  latter, 
which  integrates  thermoreceptor  output  signals  from  various  body  sites  and 
determines  the  appropriate  thermoregulatory  response,  will  not  be  discussed 
(interested  readers  should  refer  to  Stolwijk  and  Hardy  1968;  Montgomery  1974; 
Hwang  and  Rons  1977;  Arkin  and  Shitzer  1984;  Wissler  1985).  The  controlled 
system  considers  the  human  body  as  a  passive  heat  transfer  system  composed  of 
segments  representing  the  various  components  of  the  body  (e.g.,  head,  trunk,  etc.). 
Each  segment  is  further  divided  into  compartments  representing  different  tissue 
mass.  Each  compartment  is  usually  assigned  a  unique  set  of  physiological, 
thermophysical  and  geometric  characteristics  that  determine  its  temperature  through 
the  bioheat  equation  (see  Equation  1).  Figure  1  shows  the  physical  representation 
of  the  human  body  modelled  by  Stolwijk  (1970).  This  model  and  its  modification 
for  cold  water  immersion  (Montgomery  1974)  were  the  basis  for  the  present  study. 
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Pennes  (1048)  outlined  the  bioheat  equation  to  predict  the  steady-state 
temperature  distribution  of  a  compartment.  In  cold  water  immersion,  most  studies 
are  concerned  with  the  transient  or  unsteady  temperature  distribution.  Analytic 
solutions  of  the  unsteady  temperature  distribution  (Shitzer  1085;  Klinger  1085)  are 
unwieldy  for  thermoregulatory  modelling.  To  simplify  the  solution  of  the  unsteady 

I 

problem,  finite  difference  methods  are  often  applied.  Although  analyses  have  been 
conducted  for  specific  conditions  using  this  approach,  no  general  optimization  has 
been  undertaken.  The  simpler,  but  related,  heat  conduction  problem  can  be  used 
to  outline  an  approximate  optimization  scheme  for  using  the  finite-difference  method 
of  solution. 

Initial  conditions  usually  assume  a  deep  core  temperature  (T£)  and  a  metabolic 
rate  (MR)  based  on  "standard”  man  in  a  state  of  thermal  neutrality  (i.e.,  zero  heat 
storage).  This  makes  subsequent  comparison  between  predicted  and  measured  values 
of  these  variables  during  immersion  ambiguous  since  a  person's  pre-immersion  values 
do  not  necessarily  match  "standard”  man's,  especially  if  the  person  is  lean  or  fat. 
A  procedure  to  match  the  initial  values  of  T  and  MR  under  a  condition  of 
thermal  neutrality  is  outlined. 

Conventional  methods  of  predicting  the  convective  heat  loss  are  sensitive  to  the 
assumptions  of  body  shape  and  water  motion,  and  especially  to  the  skin-water 
temperature  difference  when  the  skin  temperature  (T^)  is  near  or  at  its  steady- 
state  value.  An  alternative  approach  which  relies  on  the  heat  balance  of  the  skin 
compartment  is  developed. 

Finite-Difference  Method 

The  bbheat  equation  expressed  in  cylindrical  co-ordinates  is 


(pcV) 
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where  r  is  the  radial  distance  from  the  center,  k  is  the  thermal  conductivity,  T  is 
th  '  iperature,  p  is  the  density,  c  is  the  specific  heat,  m  is  the  metabolic  rate 
per  unit  volume,  t  is  the  time,  and  the  subscript  bl  refers  to  the  blood. 
Refinements  such  as  separate  heat  exchange  with  the  arterial  and  venous  blood, 
tangential  and  angular  temperature  gradients,  and  counter-current  heat  exchang: 
(Kuznetz  1079;  Arkin  and  Shitzer  1984;  Wissler  1985)  are  not  considered  here.  To 
test  the  accuracy  of  the  finite-difference  solution  of  Equation  1,  predictions  are 
either  compared  to  experimentally-determined  values  (Chan  et  al.  1973)  or  to  exact 
solutions  that  are  made  possible  by  the  selection  of  specific  initial  and  boundary 
conditions.  One  such  set  to  predict  temperatures  in  living  tissue  was  used  by 
Gordon  and  Roemer  (1975)  and  they  found  that  errors  were  reduced  by  decreasing 
the  integration  time  interval  from  10  to  1  min,  by  increasing  the  number  of  nodes 
from  4  to  10,  and  by  decreasing  the  nodal  spacing  for  large  temperature  differences. 
For  tissue  characteristics  and  conditions  different  from  the  ones  considered  by 
Gordon  and  Roemer,  a  more  general  set  of  guidelines  is  desirable  to  optimize  the 
numerical  integration  of  Equation  1. 

An  approximate  optimization  scheme  can  be  outlined  by  considering  the 
simpler  heat  conduction  problem  with  tissue-blood  heat  conduction  and  tissue 
metabolism  removed  so  that  an  exact  solution  of  the  average  tissue  temperature 
(Tavg)  c%a  k*  found.  Tav^  serves  as  a  good  reference  point  for  comparisons 
between  the  exact  and  finite-difference  solutions  of  Equation  1  since  the  change  in 
TftVg  is  directly  related  to  the  tissue  hen  storage  which  is  the  quantity  that 
determines  net  heat  exchange. 

Consider  the  problem  defined  in  Figure  2  where  the  initial  tissue  temperature 
has  a  uniform  value  a  and  the  surface  temperature  (T#)  is  fixed.  The  exact 
solution  of  T4Vg,  normalised  for  dimensionless  analysis,  is  given  by  (Carslaw  and 
Jaeger  1980;  Crank  1975) 


4 


T„,  -  I.  -  1  -  £^-i-)2'  *xp  C'  K“iSt) 


where  a.  is  the  positive  root  of  Jo  (a.R),  Jo  is  the  Bessel  function  of  the  first  kind 

of  zero  order,  end  K  equals  k/ pc.  Values  of  Equation  2  are  plotted  against  the 
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dimensionless  quantity  Kt/R  in  Figure  3. 

The  explicit  finite-difference  solution  of  Equation  1  using  a  three-point 
approximation  is  given  by  (Ozisik  1985);  for  the  innermost  nodal  point  (n=0) 

T  (t  ♦  At)  *  T  (t)  +  JCAt*  (Tj  (t)  -  T  (t))  (3a) 


and  lor  1  <=  n  <=  N-l 


T  (t  ♦  At)  *  T  (t) 


/W  •  <T: 
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\Vis  %y  _ 
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where  At  is  the  integration  time  interval,  r  and  q  are  radii  defined  in  Figure  2, 
and  the  subscript  refers  to  the  nodal  point.  The  outermost  nodal  point  (n=N) 
temperature  is  fixed  at  T  .  T  is  determined  by  weighting  the  temperature  of 
each  node  according  to  the  thermal  capacity  of  tissue  around  the  node  defined  by 
the  boundaries  at  q.  Values  of  normalized  TfcV^  for  3,4,6  and  10  codes  spaced 
apart  by  equal  volumes  are  shown  in  Figure  3.  These  values  were  obtained  by 
choosing  At  cs  IQ**  R‘/K  which  allowed  ten  integration  intervals  prior  to  the 
earliest  values  shown. 

Four  factors  contribute  to  the  size  of  error  using  the  finite-difference  solution: 
they  are  l)  value  of  At;  2)  number  of  nodes;  3)  time  of  prediction;  and  4)  nodal 
spacing.  Decreasing  At  further  than  its  present  value  led  to  improvements  of  less 
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than  0.1%.  Therefore,  the  size  of  error  indicated  in  Figure  3  is  only  minimally 
attributed  to  the  present  >ce  of  At.  Increasing  At  worsens  the  prediction  but  the 
extent  depends  upon  the  number  of  nodes  used.  As  shown  in  Table  1,  the 
prediction  is  increasingly  sensitive  to  At  as  the  number  of  nodes  increase. 

An  infinite  number  of  nodes  will  model  the  heat  conduction  problem  exactly, 
therefore,  increasing  the  number  will  logically  decrease  the  error  as  shown  in  Figure 
3.  The  minimum  number  of  nodes  required  by  the  present  finite-difference  method 
is  three,  yet  after  a  long  period  of  time  the  error  using  three  nodes  becomes 
negligible.  Therefore,  the  optimal  choice  of  the  number  of  nodes  depends  upon  the 
acceptable  size  of  error  at  a  given  time.  This  can  be  best  illustrated  by  a  practical 
example.  Consider  a  cylindrical  model  compartments  of  radius  R  =  2  inch  (5  cm) 
and  a  typical  tissue  conductivity  K  =  0.00016  inch  /s  (0.001  cm  /s).  If  the  error 
must  not  exceed  10%  after  1  min  (i.e.,  at  Kt/R  =  0.0024),  then  ten  nodes  can  be 
used  as  seen  in  Figure  3.  If,  however,  this  constraint  is  relaxed  so  that  the  same 
level  of  error  must  not  be  exceeded  after  10  min,  then  four  nodes  can  be  used.  Of 
course,  this  example  is  predicated  on  the  use  of  At  =  10  R  /K,  hence  some 
adjustment  may  be  necessary  if  At  is  increased.  Otherwise  Figure  3  is  a  reliable 
guide  to  approximate  the  optimal  number  of  nodes  for  any  combination  of  tissue 
conductivity  and  cylindrical  dimension. 

The  size  of  error  using  a  geometry  that  separates  nodes  by  equal  volumes 
(EV),  as  in  the  present  study,  is  considerably  smaller  than  that  generated  by  using 
equally  spaced  nodes  (EN).  This  is  demonstrated  in  Table  1  where  the  percent 
error  is  listed  for  both  ways  of  separating  nodes.  This  result  is  consistent  with  the 
finding  of  Gordon  and  Roemer  (1975)  that  the  error  is  reduced  by  decreasing  the 
nodal  spacing  where  temperature  gradients  are  highest. 
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These  guidelines  are  generally  applicable  for  the  heat  conduction  problem 
outlined  in  Figure  2.  The  extent  w  ich  these  guidelines  would  change  with  the 
inclusion  of  tissue-blood  heat  conduction  and  tissue  metabolism  is  not  known.  We 
suggest,  therefore,  that  the  above  guidelines  serve  as  an  approximate  starting  point 
for  optimizing  the  finite-difference  solution  of  the  unsteady  bioheat  equation. 

Initial  T  and  MR  values 

. . — '‘C"  -  '  ■  . . 

Initial  conditions  usually  assume  a  state  of  thermal  neutrality  which  is  the 
starting  point  for  many  experimental  studies.  The  conventional  procedure  for 
determining  the  model’s  initial  temperature  distribution  is  to  assume  the 
physiological „  thermophysical  and  geometrical  characteristics  of  "standard'*  man  and 
simulate  an  exposure  to  air  within  the  zone  of  thermal  neutrality  (Stolwyk  1070). 
The  thermoregulatory  model  will  then  generate  a  steady-state  temperature 
distribution  with  zero  heat  storage  in  each  compartment.  However,  when  comparing 
model  predictions  to  actual  measurements,  the  initial  values  of  T£  and  MR  between 
the  model  and  the  subject  may  not  be  matched.  Unless  these  values  are  matched, 
it  would  be  difficult  to  ascertain  to  what  extent  differences  between  these  values 
during  the  subsequent  immersion  period  are  attributed  to  the  model.  In  addition, 
initial  differences  in  Tc  and  MR  between  the  model  and  the  subject  impose  an 
arbitrary  thermoregulatory  response  at  the  outset  of  the  exposure.  Simply  applying 
a  subject's  anthropometric  characteristics  to  the  model  using  the  conventional 
procedure  will  not  assure  that  initial  values  of  ?c  and  MR  between  model  and 
subject  are  matched.  The  following  describes  a  procedure  for  matching  these  values 
under  a  condition  of  thermal  neutrality,  similar  to  the  method  used  by  Miller  and 
Seagrave  (1074)  in  another  application. 

Consider  a  model  consisting  of  i  number  of  segments  each  consisting  of  J 
number  of  compartments.  The  heat  storage  of  the  jth  compartment  of  the  ith 
segment  is  given  by  Stolwjjk  (1970). 
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HFiJ  =  QiJ  '  Ei j  '  Bcij  +  TDiJ-l  "  TDiJ  W 

where  Q  is  the  metabolic  rate,  E  is  the  ev  '  tive  heat  loss,  BC  is  the  heat 

conducted  to  the  central  blood,  and  TD  is  the  heat  conducted  to  the  adjacent  outer 

compartment.  The  latter  two  terms  are  explicitly  given  by, 

BCy  =  VBIy-(TM  •  Tbi>  W 

and 


TD. .  =  TC.  .*(T. .  -  T. . . ,)  (6) 

ij  ij  1  »j  ij+1'  v  1 

where  c^j  is  the  volumetric  heat  capacity  of  the  blood,  BF  is  the  blood  flow  rate, 
and  TC  is  the  thermal  conductance  (defined  by  the  coefficients  in  Equation  3). 
For  the  innermost  compartment. 


HFi.l  =  <\l  -  Ei,l  ■  DCU  -  TDi.i  <7> 

where  E  represents  the  respiratory  heat  loss  for  the  head  and  trunk  segments  only, 

otherwise  its  value  is  zero.  For  the  outermost  compartments, 

HFu  =  •  Ei.J  -  BCij  +  TOu-i  •  Hi*(Tu  •  *«*>  <e> 

where  H  is  the  combined  radiative  and  convective  heat  transfer  coefficient  in  air, 
and  T  ^  is  the  ambient  temperature.  For  the  central  blood. 


HF, 


bl 


=  t 

i-1  1-1 


(9) 


Under  a  condition  of  thermal  neutrality,  the  heat  storage  of  each  compartment, 
including  the  central  blood,  is  set  equal  to  zero.  An  iterative  procedure  is  used  to 
solve  the  thermal  neutral  temperature  distribution.  By  equating  the  heat  storage  of 
each  compartment  to  zero  and  choosing  a  starting  value  of  7-0j,  the  a>v«em  of  heat 
storage  equations  reduces  to  I  independent  sets  of  J  equations  each  that  can  be 
represented  in  matrix  notation.  For  example,  the  matrix  for  a  segment  -omprised 
of  four  compartments  (J  =  4)  as  (illustrated  in  Figure  1)  is  given  by 
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bn*i,i 

4  ICi4)  -  TCi  2  0 

*  1 

-tcm 

(^•BPi  ^  ♦  TGi  i  ♦  TCi|3)  “  TCi,3  •'  0 

0 

’  TCi,2  (cbl#BFi,3  *  TCi,2  4  ~  ^1,4 

0 

P  -  I0i,3 

*  Ki,3  *  V  J 

*T  ~ 

Ti»l 

\i  -  Ei,i  *  cw,®;y,bi 

Ti,2 

^i,2  *  2 

(10) 

Ti,3 

«i,3  *  eU’BFi.3*Tbl 

Ti,4 

^i,4  '  Bi,4  *  cbl-BFl,t'Tbl  4  Bi*Tiab 

For  all  compartments  except  that  representing  the  muscle,  Q  is  given  by  the 
compartment's  basal  metabolic  rate.  For  the  muscle  compartment, 

Q.  2  =  (IMR  -  BMR)*WORKM-  +  QB.^  (U) 

where  IMR  and  BMR  are  the  subject's  initial  and  basal  metabolic  rates  and 
WORKMj  is  the  fractional  contribution  of  the  ith  segment  to  the  difference  between 
IMR  and  BMR  (Stolwjjk  1070).  The  Gauss  Elimination  Method  (Kreyssig  1007) 
can  be  used  to  solve  for  the  temperatures  T-  T,  ,,  T,  -  and  T;  The  iteration 
continues  until  values  of  T^j  and  are  found  such  that  the  trunk  core 

temperature  of  the  model  equals  the  subject's  deep  body  temperature  (usually  the 
rectal)  and  the  heat  storage  of  the  blood  equals  sero.  The  temperature  distribution 
obtained  for  'standard*  man  by  this  procedure  agrees  to  that  obtained  using  the 
conventional  procedure,  but  for  conditions  other  than  "standard*,  this  procedure  has 


the  advantage  of  matching  model  and  subject  initial  values  of  T  and  MR. 

As  an  exampb  -jf  the  differences  that  can  be  expected  between  individuals  of 
different  morphology,  consider  the  thermal  neutral  temperature  distributions  obtained 
under  the  same  basal  am.  environmental  conditions  for  a  lean  and  a  fat  individual 


0 


shown  in  Figure  4.  Because  of  'he  higher  body  insulation  of  the  fat  individual, 
internal  temperatures  are  higher,  yet  skin  temperatures  are  close  wv  3e  of  the 
lean  individual.  In  other  words,  the  fat  individual  has  a  higher  temperature 
difference  between  the  inner  and  outer  compartments  than  the  lean  individual  untisr 
the  same  thermal  neutral  condition.  Therefore,  the  model  would  assign  different 
thermoregulatory  set-points  for  these  individuals. 

Convective  Heat  Loss 

The  conventional  prediction  of  convective  heat  loss  (C)  in  water  is  high* 
sensitive  to  the  skin-water  temperature  difference,  especially  as  Tgj{  nears  ta* 
temperature  of  the  water  (Tw): 

C  =  hc.(T  t  -  TJ  (12) 


where  h  is  the  convective  heat  transfer  coefficient.  In  nude  immersion,  the  steady- 

state  skin  temperature  (T  ,  )  is  usually  0.5  to  1.5°  C  above  T  .  Small  variations 

w 

in  T.k  about  the  value  of  Tg^  can  therefore  cause  large  variations  in  the 
predicted  value  of  C.  Such  variations  in  C  critically  affect  the  heat  storage  of  the 
akin  compartment  which  in  turn  can  perpetuate  a  fluctuation  in  the  predicted  value 
of  T^.  Unless  the  integration  time  interval  is  sufficiently  small,  these  fluctuations 
may  be  unacceptably  large  and  instability  may  result.  This  problem  is  exacerbated 
by  the  assumptions  of  body  shape  and  water  motion  that  determine  h£  (Sekias  and 


Emery  1982)  which  partly  explains  the  wide  disparity  among  values  reported  in  the 


literature  (Boutetier  et  ai.  ISTT).  To  circumvent  this  difficulty  of  a  potentially 
fluctuating  value  of  C,  an  alternative  method  of  predicting  this  value  Is  proposed, 


as  follow*. 

It  has  been  shown  experimentally  that  Tgj.  of  a  nude  subject  falls 
exponentially  during  immersion  in  cold  water  (Strong  et  al.  1984).  Therefore,  the 
rate  of  change  of  T,i  should  be  proportional  to  the  difference  between  T.i  and  its 
current  value,  expressed  in  differential  form  as 


S3 
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**  =  X*(Tsk  -Td.)  W) 

£8 

where  X  is  a  proportionality  coefficient.  The  change  in  T,s1,  over  a  time  interv~’  ' 
in  which  X  is  fixed  b  given  by 

AT*  =  (Tsk  -  T  k  )•(’.- «p(-X.A‘))  (14) 

SS  O 

where  Tg^  is  the  skin  temperature  at  the  start  of  the  time  interval.  At  the  same 
time,  the  rate  of  change  of  Tg^  must,  by  definition,  equal  the  skin  compartment’s 
heat  storage  (HFg^)  divided  by  its  heat  capacity  (Cg,^).  By  equating  this  quantity 
to  Equation  13,  the  value  of  X  is  defined  which  for  small  values  of  At  can  be 
approximated  by: 


X=  HVC6k*<Tsk  -  Tsk>  <15> 

ss  o 

During  the  unsteady  state  of  falling  Tg^,  the  value  of  HF^  is  determined 
through  Equation  8  with  the  modification  that  the  last  term  is  replaced  by  the 
exDression  for  C  given  by  Equation  12.  Although  the  sensitivity  discussed  above 
with  predicting  C  is  present  during  this  period,  this  sensitivity  diminishes  as 

falls.  Once  Tg^  nears  Tg^  (a  difference  of  0.005°  C  is  considered  negligible),  the 

as 

value  of  C  can  be  predicted  by  assuming  zero  heat  storage  of  the  skin  compartment 


(i.e.,  HF  ,  =  0)  so  that 


sk 

ci  =  <JW 


BC.  j  +  TD.  T  , 
i,J  i,J-l 


(16) 


In  this  way,  the  value  of  C  during  steady-state  of  the  skin  temperature  is  not 

sensitive  to  the  skin-water  temperature  difference. 

Throughout  this  development,  a  value  of  has  been  implicitly  assumed. 

ss 

Although  the  value  of  T  ,  cannot  be  wholly  predicted  in  advance,  a  good 

*  s  s 

estimate  can  be  made  because  of  the  narrow  range  of  values  found  experimentally 
for  nude  immersion  in  water.  In  addition,  the  values  of  BCh  j  and  TD.  are  not 
sensitive  to  the  choice  of  Tg^  ,  so  little  error  can  be  expected  using  Equation  16 
during  steady-state  of  the  skin  temperature. 
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Conclusions 

Three  specific  problems  with  modelling  the  controlled  system  of  human 
thermoregulation  for  cold  water  immersion  have  been  examined:  1)  the  finite- 
difference  solution  of  the  bioheat  equation;  2)  differences  between  predicted  and 
measured  values  of  Tc  and  ME;  and  3)  the  prediction  of  C.  An  approximate 
optimization  scheme  for  using  the  finite-difference  solution  of  the  bioheat  problem 
was  outlined  by  considering  the  simpler,  but  related,  general  heat  conduction 
problem.  In  concurrence  with  other  studies,  it  was  found  that  the  error  was 
reduced  by  decreasing  the  integration  time  interval,  by  increasing  the  number  of 
nodes,  and  by  decreasing  the  nodal  spacing  for  large  temperature  differences. 
However,  by  a  judicious  choice  of  these  factors,  the  finite-difference  procedure  can 
be  optimized.  For  example,  a  greater  benefit  is  obtained  by  doubling  the  number 
of  nodes  rather  than  halving  the  integration  time  interval. 

To  edict  a  subject’s  initial  temperature  distribution  under  a  condition  of 
thermal  neutrality,  zero  heat  storage  was  assumed  for  each  model  compartment  and 
the  refulting  system  of  linear  heat  balance  equations  were  solved  simultaneously. 
An  iterative  procedure  was  adopted  so  that  the  solution  yielded  a  value  of  trunk 
core  temperature  equal  to  the  subject’s  deep  core  (rectal)  temperature.  Since  set- 
point  values  for  thermoregulation  are  based  on  the  thermal  neutrai  temperature 
distribution,  this  procedure  of  matching  the  model  and  subject  initial  temperatures 
avoids  any  arbitrary  thermoregulation  at  the  outset  of  the  exposure. 

To  circumvent  the  uncertainties  associated  with  predicting  the  convective  heat 
loss  during  steady-state  of  the  skin  temperature,  a  method  was  developed  to  assure 
a  smooth  exponential  fall  in  Tgj.  to  a  steady-state  value.  Once  this  value  was 
attained,  zero  heat  storage  of  the  skin  compartment  was  assumed  so  that  C  could 
be  predicted  through  a  heat  balance  of  the  skin  compartment.  In  this  way,  the 


predicted  value  of  C  is  not  sensitive  to  the  small  skin-water  temperature  difference 
that  occurs  with  nude  immersion  in  water.  This  is  potentially  more  important  for 
exposures  less  severe  than  cold  water  immersion  since  the  relative  error  using 
incorrect  set-points  is  greater .  when  actual  temperatures  are  close  to  neutral 
conditions. 

Although  the  problems  addressed  in  this  paper  arose  out  of  a  need  for  greater 
accuracy  when  modelling  the  human  thermoregulatory  response  for  immersion  in  cold 
water,  the  procedures  for  handling  these  problems  may  be  applicable  in  other 
circumstances,  especially  where  finite-difference  methods  are  used  and  matching  of 
initial  conditions  are  important. 


The  views,  opinions,  and/or  findings  contained  in  .  this  report  are  those  of  the 
authors  and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,  policy,  or  decision,  unless  so  designated  by  other  official  documentation. 
Approved  for  public  release;  distribution  unlimited. 
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Table  1:  Percent  error  at  Kt/R*  =  0.01  for  nodes  spaced  apart  by  equal  volumes 
(EV)  and  by  equal  radial  displacement  (EN) 


Number  At  =  0.0001  R*  /  K  At  =  0.001  R*  /K 

of  nodes 


EV 

EN 

EV 

EN 

3 

41.0 

117.2 

41.0 

117.2 

4 

17.7 

64.2 

18.1 

64.2 

6 

6.0 

25.1 

7.0 

27.2 

10 

1.0 

7.4 

3.3 

8.4 
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FIGURE  CAPTIONS 

Figure  1.  Schematic  of  the  human  body  (one  tide  shown)  modelled  by  Stolwijk 
(1070).  Each  body  segment  is  composed  of  four  concentric  annular  compartments, 
the  head  modelled  as  a  sphere  and  the  others  as  cylinders.  Length  of  the  cylinders 
are  given  in  inches.  Radial  dimensions  are  determined  by  the  mass  of  the 
individual  modelled. 

Figure  2.  Schematic  of  the  heat  conduction  problem. 

Figure  3.  Normalized  values  of  the  average  tissue  temperature  plotted  against  the 
dimensionless  quantity  Kt/Ra.  The  solid  line  represents  the  exact  solution  of  the 
heat  conduction  problem  outlined  in  Figure  2  and  the  dashed  lines  represent  the 
finite-difference  solution.  Numbers  above  the  dashed  lines  indicate  the  number  of 
nodes  used.  The  integration  time  interval  was  R3/K  x  10'*. 

Figure  4.  Thermal  neutral  temperature  distributions  for  the  Stolwijk  (1070)  model 
assuming  basal  conditions  (BMR  =  14.5  BTU/ft*/h  (45.7  w/ma))  and  a  simulated 
exposure  to  a  86°  F  (30°  C)  environment  with  a  20  ft/min  (10  cm/s)  air  movement. 
The  open  bars  indicate  the  temperature  distribution  of  a  lean  person  (height  =  68 
in  (173  cm),  weight  =  140  lb  (C3.5  kg),  body  fat  =  7.3%)  and  the  hashed  bars 
indicate  the  temperature  distribution  of  a  fat  person  (height  =  68  in  (173  cm), 
weight  =  200  lb  (00.7  kg),  body  fat  —  23.8%).  Compartments  of  each  segment  arc 
core,  muscle,  fat  and  skin  ordered  from  left  to  right. 
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